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There is a general agreement that galaxy structures exhibit fractal properties, at 
least up to some small scale. However the presence of an eventual crossover towards 
homogenization, as well as the exact value of the fractal dimension, are still a 
matter of debate. I summarize the main points of the this discussion, considering 
also some galaxy surveys which have recently appeared. Further I discuss the 
implications for the standard picture of the observed fractal behaviour in galaxy 
distribution. In particular I consider the co-existence of fractal structures and 
the linear Hubble-law within the same scales. This fact represents a challenge 
for the standard cosmology where the linear Hubble law is a strict consequence of 
homogeneity of the expanding universe. Finally I consider the comparison of CDM- 
like models with the data noting that the simulations are not able to reproduce 
the observed properties of galaxy correlations. 



1 Introduction 

It has been well known for over twenty yearslffl that galaxy structures exhibits 
fractal properties at small scales (~ lQh~^Mpc). The scale invariant correlated 
behavior corresponds to the existence of large scale structures (hereafter LSS) 
in the galaxy distribution. This evidence came from the analysis of the angular 
galaxy catalogs and from few sparse measurements of redshifts. At scale larger 
than ^ lOh^^Mpc it was reasonable, from the angular data alone, to assume 
the homogeneity of galaxy distribution, that is the cornerstone of the Big Bang 
model. 

However, with the extensive measurements of redshifts started in the eight- 
ies it was discovered that the extent of galaxy structures and voids is limited 
only by the size of the available samples. Such a situation has been confirmed 
by several recent 3-d galaxy catalogs. It is important to note that at the present 
time the investigation of the large scale distribution of galaxies has having an 
exponential growth which will lead, in less than ten years, to collect more than 
one million of galaxy redshifts. Such a situation, together with the use of the 
modern methods of statistical physics for a quantitative characterization of the 
distribution, gives rise to a rather different perspective on the properties of the 
large scale distribution of matter in the universe as well as on the theoretical 
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methods adopted. Very recently a wide debate on this subject is in progress: 
see the web page http://www.phys.umromal.it/DOCS/PIL/pil.htmi where ah 
these materials have been collected. 



2 The problem of the usual perspective 

The usual analysis! §i of galaxy distribution identifies a correlation length of 
about ~ 5h~^Mpc. Such a length scale should characterize the distance at 

fhich the density fluctuations is of the order the average density. Pietronero 
criticized this result on the basis that such a small characteristic length is 
actually inconsistent with the existence of LSS and huge voids larger more 
than one order of magnitude in size. The problem of the standard analysis 
of galaxy correlation lies in the a priori assumption of homogeneity. In other 
terms one usually defines an average density (n) in a given galaxy sample, and 
then one compare the density fluctuations 6n to such a value. Such a procedure 
does not allow one test whether the average density is a meaningful quantity. 
However from the studies of irregular systems 13 we have learned that in a 
self similar structure there are no characteristic values and concepts like the 
average density cannot be defined properly. More specifically such quantities 
are not related to the nature of the distribution, rather they depend on the 
size of the sample, the unique meaningful length-scale which can be defined in 
such a situation. 



3 Methods of analysis 

Let us briefly illustrate the methods of statistical analysis usually used in the 
studies of irregular, self-similar, structures, but that can be successfully used 
also for the characterization of regular systems. PietroneroH proposed to study 
the conditional average density defined as 

r( ) - ^ dN _ BD 
S{r) dr Att 

Such a quantity measures the number of points dN in a spherical shell of 
thickness dr and volume S{r) = Anr'^dr, located at distance r from an occupied 
point. Then one determines the average over all the points contained in a given 
sample. Being an average quantity, T(r) is a rather stable and robust statistical 
indicator. The last equality in eq.^ holds for a fractal with dimension D < 3 
and prefactor B. If the distribution is homogeneous {D = 3) r(r) equals the 
average density in the sample. The power law behaviour of r(r) implies the self- 
similarity of the distribution. The prefactor B can be defined for real structures 
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contained in finite samples. It gives the normalization of the amplitude of the 
space density. A very simple interpretation of such a quantity is the following: 
in a homogenous distribution (according to eq.Q) B is the space density (a part 
trivial prefactor). The average distance between nearest particles is known to 
be A ~ ^1/3. In the case of a fractal distribution it is possible to showQ that 
the average distance between nearest neighbors is of the order A ~ B^/^ (this 
is an intrinsic quantity that does not depend on sample size). Note that while 
in a homogeneous distribution A gives also a reasonable order of magnitude 
for the typical voids contained in the sample, in the fractal case the size of the 
voids scales as a function of the size of the sample and it is closely related to 
another property called lacunaritp u. 

From eq.0 it follows that for a fractal distribution, the average density in 
a sample of radius Rg scales as R^~^ and hence it does not represent a mean- 
ingful reference value. This simple observation shows that the usual statistical 
analysis based on concepts like £,{r), the power spectrum and other related 
quantities becomes meaningless, unless a very well define transition to homo- 
geneity is present in the sample. In terms of r(r) this transition should be 
shown by the break of the power law into a flatter behavior with scale. 

It is simple to showLD that all the characteristic length scales usually iden- 
tified in the study of the LSS become spurious and dependent on the sample 
size, unless the density shows a constant behaviour with scale. Some example 
are: the so-called correlation length ro (scale at which ^(ro) = 1), the turnover 
scale As of the power spectrum P(fc) (defined ad dP{Xs{k))/dk — 0), the scale 
at which the density fluctuaidpns is SN/N — 1. (For a more detailed discussion 
see Sylos Labini et al, 19980). 

4 Review of main results 

A real fractal structure can be observed and defined in finite samples. Hence 
it is important to clarify the lower and the upper cut-offs among which its 
properties can be properly studied by a suitable correlation analysis. As al- 
ready mentioned the lower cut-off A is an intrinsic quantity of a given fractal 
and it gives the order of magnitude of the average distance between nearest 
neighbors. The upper cut-off Rg is defined to_he the radius of the maximum 
sphere fully contained in the sample volume El Q. This definition avoids any 
assumption on the treatments of the boundary conditions of the sample in the 
correlations analysis. 

In the range of scale A < r < i?s we have found, in a complete and 
extensive study of all the available redshift catalogs D, that galaxy distribution 
has well defined fractal properties. In particular, by the full correlation analysis 
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we have found a fractal dimension D = 2.1 ± 0.2 in the range of scales 0.5 < 
r < 150h~^Mpc (Figj^). This result is substantially stable in the different 
catalogs we have considered. A recent analysis of the SSRS2 and CfA2 South 
galaxy sampleftj is in complete agreement with all the other galaxy surveys 
previously considered. Also the ESP survejU shows a continuation of the 
fractal behavior with D « 2 up to ~ 600 800h~^Mpc, although this result is 
much weaker from a statistical point of view (see Joyce et ai, 1998 for a more 
detailed discussion on this subject). 
The prefactor B is found to be 

B = 10^15h-^Mpc-" (2) 

and hence the lower cut-off is of the order A ^ 0.5h~^Mpc. This is the min- 
imum statistical distance beyond which the statistical properties are well de- 
fined. However in real samples on should consider also a luminosity selection 
effect which can cause to change A by more than one order of magnitude, 
depending on the volume limited sample considered. 



5 Problems of the Standard FRW scenario 

In his classic paper, HubblS found a roughly linear relation between the 
spectral line displacement z = {Xobs — ^em)/^em of the line emitted by a far 
away galaxy Aem, and its distance r. The empirical Hubble Law may be written 
as 

cz = Hor , (3) 

where c is the velocity of light and Hq is the Hubble constant. As an obser- 
vationally established relation, the Hubble law does not refer to any interpre- 
tation of redshift. Space expansion and Doppler mechanisms in fait space for 
redshift, yields at first order to Eq.^. If redshift is interpreted as a motion 
effect, then 

V « Hor (4) 

where V is either space expansion velocity v^xp or ordinary velocity of a body 
moving in the Euclidean space. Usually this velocity-distance relation is called 
the Hubble Law, but it is more correct to regard it as the redshift-distance 
relation of Eq.^. This is based on the primarily measured quantities (redshift 
and distance), while velocity is inferred from redshift in the frame of some 
cosmological model. 

Since its discovery, the validity of the Hubble law has been confirmed 
in an ever increasing distance interval where local and more remote distance 
indicators may be tied together. Recently, several new distances have been 
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Figure 1: Full correlation analysis for the various available redshift surveys in the range of 
distance 0.5-^ 1050/i~^Mpc. The slope is —1, which corresponds to fractal dimension D = 2. 

(From Sylos Labini et al, 1998). 
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measured to local galaxies using observationsjif Cepheid variable stars, thanks 
to the Hubble Space Telescope programmeO. Along with previous Earth- 
based Cepheid distances, methods like Supernovae la and TuUy-Fisher have 
been better calibrated than before and confirm the linearity with good accuracy 
up to z « 0.1. Brightest cluster galaxies trace the Hubble law even deeper, 
up to Z|_ai 1, and radio galaxies have provided such evidence at still larger 
redshifto. 

It is well known that there are small deviations SV from the Hubble velocity 
Vh, connected with local mass concentrations such as the Virgo Cluster, and, 
possibly the Great Attractor. However, these perturbations are still only of 
the order SV/Vh ~ 0.1, while in the general field the_pubble law has been 
suggested to be quite smooth, with SV around 50km /xl. 

Without the actual knowledge of matter distribution, the linearity and the 
small scatter of the observed Hubble law for field galaxies would make one guess 
that the galaxies are uniformly distributed: as it was asserted above, this is the 
basis for the linear Hubble law in the standard cosmology. In fact, it has been a 
common supposition that when the Hubble law was found in the nearby space, 
one finally had entered a cosmologically representative region of the Universe. 
At the same time, it has been clear that at small distances where Hubble found 
his relation, the galaxy distribution is quite inhomogeneous. Though, it has 
been believed that beyond some, not too large distance, the distribution should 
become uniform. 

As we have already mentioned, studies of the 3-dimensional galaxy uni- 
verse have shown that de Vaucouleurs' prescient view on the matter distri- 
bution is valid at least in the range of scales ^ 0.5 200 Mpc (hereafter 
Hq = 55km/ sec/ Mpc). The Hubble and de Vaucouleurs laws describe very 
different aspects of the Universe, but both have in common universality and 
observer independence. This makes them fundamental cosmological laws and 
it is important to investigate the consequences of their coexistence at similar 
length-scales. In Fig|| we display these laws together. A_representative Hub- 
ble law has been taken from Fig. 4 of Teerikorpi (1997^, based on Cepheid 
distances to local galaxies, Tully-Fisher distances from the KLUN programme, 
and Supernovae la distances. The behavior of the conditional density (De Vau- 
couleurs law) presented in Fig.|^ has been taken from Sylos Labini et aZ.(1998). 

The puzzling conclusion from Fig.|^ is that the strictly linear redshift- 
distance relation is observed deep inside the fractal structure. (Note that 
in the analysis of galaxy redshift surveys one uses the Hubble law for the 
distance determinations as an experimental fact, i.e. any assumption has been 
used) . This empirical fact presents a profound challenge to the standard model 
where the homogeneity is the basic explanation of the Hubble law, and "the 
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Figure 2: Hubble redshift-distancc and dc Vaucouleurs density-distance laws in the distance 
scales from 1 to 500Mpc {Hq = bbkm/ s/Mpc). The Hubble law (increasing from left to 
right) is constructed from: galaxies with Cepheid-distances for cz > (triangles), galaxies 
with TuUy-Fisher (B-magnitude) distances (stars), galaxies with SNIa-distances for cz > 
3000fcm/s (filled circles) . TF-distance points are generally averages of a few tens of galaxies 
from the "unbiased plateau" of the method of normalized distances. Redshift cz is reduced 
to the Local Group center and contains the small correction due to the Virgo infall velocity 
field. The solid line corresponds to the Hubble law with Ho = bbkm/s/Mpc. The de 
Vaucouleurs law (decreasing from left to right) in the normalized form is constructed from 
the computation of the conditional average density in the following redshift surveys: CfAl 
(crosses), Perseus-Pisces, LCRS (filled diamonds), ESP (triangles left) and LEDA (Sylos 
Labini et al, 1998). The normalization between the different densities takes into account 
the different magnitude limited of the various redshift surveys. The dotted line corresponds 
to the de Vaucouleurs law with correlation exponent 7 = 1, i.e. D = 2. (From Baryshev et 

al, 1998) 
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connection between honiogcncity and Hubble's law was the first success of tlie 
expanding world model" EZl. This also reminds us the natural reaction of several 
authors: "In fact, we would not expect any neat relatiDn of proportionality 
between velocity and distance [for such close galaxies]" Ej. 

However, contrary to the expectations, modern data show a good linear 
Hubble law even for nearby galaxies. How unexpected this actually is, can be 
expressed quantitatively for the standard model and is briefly discussed below 
(for a more detailed discussion see Baryshev et ai, 1998). 

According to the standard Big Bang model the universe obeys Einstein's 



Cosmological Principle: it is homogeneous, isotropic and expanding li^. Ho- 
mogeneity of matter distribution is the central hypothesis of the standard 
cosmology because it allows one to introduce the space of uniform curvature in 
the form of the Robertson- Walker line element R{t) . This line element leads 
immediately to a linear relation between velocity and proper distance. Indeed, 
consider a comoving body at a fixed coordinate distance from a comoving ob- 
server. At cosmic t, let I — R{t)y be the proper distance from the observer. 
The expansion velocity v^xp — dl/dt, defined as the rate of change of the proper 
distance /, is 

Vexp ^ HI ^ c ■ — (5) 
Ih 

where H = R/R is the Hubble constant and Ih = c/H is the Hubble distance. 
In this way, the linear velocity-distance relation of Eq.|| is an exact formula 
for all Friedmann models and a rigorous consequence of spatial homogeneity. 
In particular, for I > Ih, the expansion velocity Vexp > c. Such-axi apparent 
violation of special relativity is consistent with general relativity c2l. 

In the expanding space the wavelength of an emitted photon is progres- 
sively stretched, so that the observed redshift z is given by Lemaitre's redshift 
law 

^ Rjtobs) _ 

~ R{t 

em } 

which is a consequence of the radial null-geodesic of the FRW line element. 
For z Eq.^ yields z ~ dR/R w Hodt « 1/Ih, and from Eq.|| one gets the 
approximate velocity-redshift relation that is valid for small redshifts 

Vexp ~ CZ . (7) 

We note that the expansion velocity-redshift relation differs from the rela- 
tivistic Doppler effect. So, the space expansion redshift mechanism in the 
standard model is quite distinct from the usual Doppler mechanism. We stress 
this points, because in the literature these two redshift mechanisms are often 
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confused. In the context of the standard cosmology, it has been natural to 
interpret the Hubble Law as a reflection of Eq.|^ and to regard the coefficient 
of proportionality Hq in Eq.|| as the present value of the theoretical Hubble 
constant H from Eq.||. 

We consider now the case of an expanding universe, where an average 
density is well defined and has a constant value po. In such a case, by neglecting 
relativistic effects and the terms depending on pressure, according to the linear 
approximation, there is a velocity deflection SV from the unperturbed Hubble 
flow Vh — HoT in the scale where the density perturbation is 5p. In the case of 
zero cosmological constant and spherical mass distribution, this deflection has 
grown during the Hubble time to the present value which is (Eq.20.55 from 
Peebles, 1993): 

^ = ln,°-'^-P (8) 

Vh 3 Po 

where ilo — Po/Pcrit is the density parameter of the Friedmann model. This 
approximation holds in the limit Sp/ po <C 1. 

Let us consider a two-component model for the density distribution in 
a Friedmann universe. First, there is the component which exhibits fractal 
behavior up to a maximum scale, and which we call Aq. At larger scales 
this component is homogeneous with an average density pium- The second 
component is dark matter, homogenous at all scales, with density pdark (see 
flgj^). For such a model there is a definite constant density at all scales larger 
than Aq. This density po is the sum of pium and pdark- This means that 
the behavior of this model at scales larger than Aq is identical to that of the 
Friedmann model for which 

— ^dark ^lum ■ (9) 

At such large scales the Hubble law is unperturbed. The density distribution 
of luminous matter for scales r < Aq, can be written as 

Plum{r) = Plum (—] (10) 



r 

where 7 = 3 — D as usual. For the scales r > Aq we have that 

Plurn{r) = Plum = P/«m(Ao) . (U) 

The density contrast can be written as 

5p_ _ Plumjr) + Pdark - PO ^^^^ 

Po Po 
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Figure 3: A simple two-component model: galaxies have a fractal distribution corresponding 
to correlated density fluctuations superimposed to a uniform background. In the figure it is 
shown the behaviour of the number density of galaxies and the matter density (from Durrer 

& Sylos Labini 1998) 



In terms of the Friedmann density parameters as defined above, this becomes, 
at scales smaller than Aq: 

At scales larger than Aq the density contrast clearly vanishes. 

By using the linear approximation (Eq.^), we may obtain a rough esti- 
mation of the expected deflection from the Hubble law in the two component 
model, for the scale in which the density contrast is less than 1 . Although the 
linear approximation is valid only for 5p/pQ ^ 1, the obtained results give a 
first quantitative indication of the effects of self-similar fluctuations. Moreover 
the assumption of spherical mass distribution is a rough one, and it holds only 
for average quantities. In the case of real fractals deviation from spherical 
symmetry can play an important role, at least at small scale. Under these 
approximations, the radial velocity measured by an average observer at scale 
r < Ao is (from Eq.^ 
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Actually, this is the prediction averaged over many observers in different frac- 
tal structure points (galaxies). For any particular observer, there will be a 
deflection from this average law. 

We take the maximum scale Aq of fractality and the fractal dimension D 
from the observed cosmological de Vaucouleurs Law and calculate the expected 
deflections from the Hubble law in our two-components Friedmann model. In 
Fig.^ we show three theoretical predictions for the velocity deflection in the 
case where the observed fractal structure contains all the matter, i.e. when 
^lum = ^0- We have fixed Ao = 200Mpc and fractal dimension D = 2. In 
this case the linear approximation holds for r > 130Mpc. At smaller scales we 
should consider non-linear effects which are not simple to be treated. However 
we should expect even stronger velocity perturbations, due to the highly inho- 
mogeneous structures distribution. The predictions correspond to three values 
of the cosmological density parameter fig = 1, 0.1, 0.01. From FigJJ it follows 
that such Friedmann models, purely fractal within 200Mpc, are excluded if 
f^o ?J 0.01. This confirms the previous suggestions that small fig is needed for 
hierarchic models (Sandage, et aZ.1972 - hereafter STH). 

For instance, there is one possible way to save the Friedmann universe with 
the critical density parameter r^o = 1. It was implied already by STH that 
dark matter, uniformly filling the whole universe and decreasing the relative 
density fluctuations, could reconcile the observed fractal structure with the 
linear Hubble law. However, they did not give a quantitative estimate of the 
amount of dark matter needed. With the new data on the Hubble and de 
Vaucouleurs laws, we can derive the lower limit for the amount of the needed 
uniform dark matter. In Eq.^ we fix fio = 1 find let flium have different 
values. Fig.^ gives the developed version of the STH test, now showing that 
^dark should be larger than 0.99. If the actual maximum scale of fractality 
is larger than 200Mpc (with D = 2), then the amount of luminous matter 
may be in confiict with the Big Bang nucleosynthesis prediction for baryonic 
matter. For example, if Aq > lOOOMpc (as suggested by Sylos Labini et ai, 
1998) then flium will be probably less than 0.001. 

It should be emphasized that this estimate of the amount of dark matter 
is independent on the physics of the early universe. It also does not depend 
on the determination of mass-to-luminosity ratio of galaxies. 

6 Inconsistency of CDM models 

It is simple to see that a fractal behavior of galaxy distribution with dimension 
Z? « 2 up to, at least, ~ 50h~^AIpc is not compatible with standard CDM 
models. In fig.0 we show the behavior of the fractal dimension versus distance 



11 



10' 



10^ 



VH=H„r 
V„.s ("„=o.i) 




100 



126 
r (Mpc) 



158 



200 



Figure 4: We show three theoretical predictions for the velocity deflection in the case where 
the observed fractal structure contains all the matter, i.e. when Qium = ^0- We have fixed 
Ao = 200Mpc and fractal dimension D = 2. The predictions correspond to three values of 
the cosmological density parameter Hq = 1,0.1,0.01. (Prom Baryshev et al, 1998) 
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Figure 5: In this case the total density is Qq = 1. We show the behaviour of Vobs derived 
by the Unear perturbation approximation in the case Sp/po <S 1, for various values of flium- 
The fractal dimension of luminous matter is D = up to Aq = 200A/pc. (From Baryshev et 

al, 1998) 

in three Cold Dark Matter models of power spectra with shape and normalized 
parameters (from Wu et al., 1998). We may see that a fractal dimension of 
_D « 2 at '--^ 40 ^ bOh^^Mpc is incompatible with all the models. Probably 
by varying the parameters of the simulation (or the mixture of Hot and Cold 
Dark Matter) one may hope to obtain a better agreement. Any new survey has 
required a new adjustment of the parameters and this alone shows the internal 
problems of the standard models of galaxy formation. 

We belive that the most important theoretical consequence of our results 
is that one may shift the attention of the study from correlation amplitudes 
to correlation exponents (Sylos Labini et aZ.1998, and Durrer & Sylos Labini, 
1998). 

7 Discussion 

Investigation of the large scale distribution of galaxies in the universe is now 
in a new phase, which is characterized by new observational data and new 
methods of analysis. It has become an especially hot and debated topic in 
cosmology, because the revealed fractality contradicts Cosmological Principle 
in the sense of homogeneity but not in the sense of the equivalence of all the 
observers^. Above we have discussed the fractality and its implications for 
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Figure 6: The fractal dimension versus distance in three Cold Dark Matter models of power 
spectra with shape and normalized parameters as reported in the labels (solid lines) (from 
Wu et ai, 1998). We show also the different experimental determinations of the fractal 
dimension we have obtained. No agreement can be found at any scale. 
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cosmology. Our main conclusions are: 



Observations show that there is a fractal distribution of galaxies, having 
fractal dimension 13 w 2 in the scale range from lMpch~^Mpc to, at 
least, lOOh-'^Mpc. While there is a general agreement on the small scale 
fractal properties of galaxy distribution, the actual value of LLajad the 
eventual prese nce of an upper cut-off, are still matter of debatcElEj . (See 
the web page http :/ / www.phys . uniromal At/ DO CS/ PIL /pil.htm\ where 



all these materials have been collected). 

The traditional statistical analysis based on the assumption of homogene- 
ity (i.e. S,{r)), should be replaced by the more general methods of modern 
statistical physics. Such methods are able to characterize scale-invariant 
distributions as well as regular ones. 

An isotropic fractal distribution is fully compatible with the reasonable 
requirement of the equivalence of all the observers. Hence the Standard 
Cosmological Principle, which requires isotropy and homogeneity, may 
be replaced by the Conditional Cosmological Principle. In such a case 
the condition of local isotropy around any structure point, without the 
assumption of analyticity of matter distribution, does not imply the ho- 
mogeneity of matter distribution. 

The paradox of linear Hubble law within the fractal de Vaucouleurs 
density-distance law is sharpened with the new data: strong deflections 
from the Hubble flow are expected in the framework of the standard 
Friedmann model. 

From a developed version of the old Sandage-Tammann-Hardy test we 
derive the minimum amount of the uniform dark matter, Qdark — 0.99, 
which is consistent with the presently known Hubble and de Vaucouleurs 
laws. This result is independent of the early universe physics. If the 
maximum scale of fractality is larger than 100h~^Mpc, this test may be 
regarded as crucial for the standard cosmology. 

We have shown that a fractal behavior of galaxy distribution with dimen- 
sion 13 « 2 up to, at least, ~ 50h~^Mpc is not compatible with standard 
CDM models. 

• One the most important theoretical consequence of our results is that 
one may shift the attention of the study from correlation amplitudes 
to correlation exponents. The revision of the concept of bias (Durrer 
& Sylos Labini, 1998) is an example of such a situation. De Vega et 
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have proposed a field theory approach to the fractal structure 
of the universe. In such a model the dominant dynamical mechanism 
responsible for the scale invariant distribution is self-gravity itself. This 
model represents an interesting approach on the lines of modern statis- 
tical physics previously: it focuses on the expoenents rather than to the 
amplitudes of correlation. In the near future we are planning to study 
some experimental consequences of such a theoretical framework. 
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